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ABSTRACT 


Information  matrices  are  derived  for  estimates  of  the  range  parameters  of 
moving  targets  as  obtained  by  combining  a  priori  information  (if  available)  with 
reflected  radar  signals  observed  in  the  presence  of  additive  white  Gaussian  noise. 

The  inverse  of  the  information  matrix  provides  a  lower  bound  on  the  covariance  matrix 
of  any  unbiased  parameter  estimates.  This  bound  can  be  approached  with  a  high 
signal -to-noise  ratio  and  optimum  data  processing  (matched  filters). 

Arbitrary  frequency  modulation,  amplitude  modulation  and  target  motion  as 
well  as  various  assumptions  on  processing  the  R-F  phase  are  considered.  The 
multiple  target  case  enables  investigation  of  a  signal's  resolution  ability  as  well  as  its 
accuracy  potentials.  Results  for  a  carrier  frequency  much  greater  than  the  effective 
signal  bandwidth  are  obtained  as  a  special  case. 

A  main  purpose  of  the  report  is  the  reduction  of  the  original  radar  problem 
to  a  linear  model  which  is  equivalent  in  the  sense  of  having  the  same  information 
matrix.  These  models  provide  valuable  insight  int  o  the  relative  effects  of  multiple 
targets,  choice  of  modulation,  a  priori  information,  and  assumptions  regarding  R-F 
phase  and  bandwidth.  The  linear  equivalent  model  also  leads  to  a  valuable  computational 
algorithm  for  investigations  using  digital  or  hybrid  computers.  The  various  special 
cases  of  interest  are  obtained  by  simple  modifications  of  the  general  case  and  thus  the 
algorithm  can  provide  a  very  versatile  tool  for  evaluating  and  designing  radar  signals. 
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1.  INTRODUCTION 


We  derive  general  formulae  for  the  accuracies  with  which  the  range  parameters 
of  moving  targets  can  be  estimated  from  reflected  radar  signals  observed  in  the  presence 
of  white  noise.  The  single  target  case  has  already  been  investigated  by  many  authors. 
(See  for  example,  Refs.  1  and  2.  Reference  l  contains  an  extensive  bibliogranhy.)  In 
this  sense,  our  work  on  the  single  target  case  is  a  rehash  of  the  previous  studies  in  a 
more  general  form.  However,  our  motivation  is  not  to  write  general  formulae,  but 
rather  to  convey  a  simplified  picture  of  a  complex  problem.  We  pursue  this  goal  by 
developing  linear  parameter  estimation  models  which  are  equivalent  to  the  original 
radar  problem.  These  simple  models  delineate  the  relative  effects  of  amplitude  and 
frequency  modulations,  assumptions  on  how  the  R-F  phase  is  processed,  and  aprtort 
information.  They  provide  a  convenient  starting  point  for  the  design  of  optimum 
modulations.  Of  equal  if  not  greater  importance,  they  suggest  a  simple  computational 
procedure  for  evaluating  performance  and  hence  for  conducting  an  experimental  signal 
design. 

The  formulae  for  the  multiple  target  case  enable  investigation  of  the  resolution 
ability  of  a  radar  signal  by  determining  how  well  the  parameters  of  several  targets 
can  be  simultaneously  estimated  from  observing  the  sum  of  the  radar  signals  reflected 
from  the  various  targets.  As  in  the  single  target  case,  the  goal  is  an  equivalent  linear 
model  which  leads  to  both  understanding  and  a  convenient  computational  algorithm. 

Our  analysis  is  a  straightforward  application  of  certain  results  of  classical 
statistics;  the  Cramer-Rao  inequality  and  linear  regression  theory.  We  treat  unknown 
parameters  corresponding  to  signal  amplitude  and  phase  the  same  way  as  the  unknown 
parameters  of  the  target’s  motion. 

In  Sec.  2  we  review  the  necessary  statistical  theory  and  in  Sec.  3,  apply  it  to  a 
general  radar  model.  In  Sec.  4,  the  effects  of  unknown  phase  and  amplitude  parameters 
are  discussed  in  terms  of  a  matrix  partitioning  which  removes  these  nuisance 
parameters  from  the  problem.  Equivalent  linear  models  for  various  single  target  cases 
are  presented  in  Sec.  5  while  the  multiple  target  case  is  discussed  in  Sec.  6.  A  general 
discussion  is  given  in  Sec.  7. 
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2.  ERROR  ANALYSIS  FOR  NONLINEAR  PARAMETER  ESTIMATION 


Consider  the  following  problem.  We  observe 

z(t)  =  y(t,t°)  +  ij(t)  OstsT  (2.1) 

where  y(t,£)  is  a  known  function  of  time  t  and  the  parameter  column  vector  £°  is 

a  fixed  but  unknown  value  of  £.  rj(t)  is  a  zero  mean,  stationary,  white  Gaussian 

2  o 

stochastic  process  with  two-sided  noise  power,  a  .  Unbiased  estimates  of  £  are  to 
be  made  using  the  observations  and  whatever  a  priori  information  is  available.  We 
are  interested  in  the  covariance  matrix  of  the  errors  in  the  estimates. 

The  following  notational  convention  is  used  repeatedly.  If  x(t ,$)  is  some 
function  of  time  t  and  an  unknown  parameter  vector  then 


(2.2) 


Thus  3(x,  j£)  is  the  column  vector  formed  from  the  partial  derivatives  of  x(t,  j£)  with 
respect  to  ^  ^ ^  *  *  • »  the  elements  of  The  partial  derivatives  are  evaluated  at  some 

chosen  value  of  £  denoted  by  j£°. 

Consider  first  a  discrete  time  version  of  Eq.  (2. 1)  wherein  we  make  N 
observations 

z(n)  =  y(n,  f°)  +tj( n)  n  =  l,---,N 
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where  rj(n)  is  some  discrete  time  stochastic  process.  Let 


F(z(l) . z(N),C°)  =  F(z,£°) 


be  the  cumulative  distribution  function  of  the  observations.  Let  £  be  any  unbiased 

o  A 

estimate  of  £  ,  made  from  the  z(n),  n  =  1,  •  •  • ,  N.  Define  Z  as  the  covariance  matrix 

o  A 

of  £  -£  ;  that  is 

*  a  A  n  A 

2  =  E[  (£-£)(£-£)'] 


where  the  prime  denotes  transpose.  If  we  define 


B  = 


I 


B(log  F/£)^(log  F/£)dF(z,£) 


(2.3) 


then  under  sufficient  regularity  conditions,  the  Cramer-Rao  or  Information  inequality 
states 


Z  >  B 


-1 


(2.4) 


*  -1 

where  the  matrix  inequality  implies  Z  *  B  is  a  positive  definite  matrix.  References  3, 

4  and  5  are  three  among  many  which  derive  this  result.  (In  Refs.  4  and  5  see  discussions 

on  efficiency  of  estimation.)  When  J7(n)  is  a  zero  mean,  discrete,  stationary  white 

2 

Gaussian  process  with  variance  cr  ,  Eq.  (2. 3)  becomes 


N 


B  =  4  ^  3(y/£)3'(y/£>  • 


n=l 


(2.5) 


The  corresponding  result  for  the  continuous  time  case  of  Eq.  (2. 1)  is 


T 

B  =  —  j*  3(y/£)3#(y/£)dt. 


(2.6) 


0 
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Now  consider  an  alternate  approach  to  investigating  the  errors  associated  with 
nonlinear  parameter  estimation  problems.  A  Taylor  series  expansion  of  y(t,  j;)  in 
terms  of  the  elements  of  £  about  £°  gives, 


y(t,  £)  =  y(t,  £°)  +  ?f(y/£)A£  +  Remainder  (2. 7) 


where 


A£  =  £-£°. 
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Assume  we  can  make  an  estimate,  1,  from  our  observations  and  aprtort  information 
such  that  £  -  £  is  small  enough  to  make  the  remainder  in  Eq.  (2.  7)  negligible.  Then 
the  error  analysis  of  our  original  nonlinear  problem  is  equivalent  to  an  error  analysis 
of  the  linear  parameter  estimation  model 


z(t)  =  r'(y/£)A£  +  r?(t) 


(2.8) 


where  3(y/£)  is  a  known  function  of  time  and  A£  is  to  be  estimated.  If  maximum 
likelihood  or,  in  this  case,  unbiased  least  squares  estimates  of  A£  are  made,  the 
resulting  error  covariance  matrix,  2,  is  given  by  the  well-known  result 

2  =  B~*  (2.9) 


(For  the  discrete  time  case,  see  discussions  on  linear  regression  and  the  linear 
hypothesis  in  Ref.  3  through  6.) 

The  matrix  B  is  called  the  information  matrix*  (or  Fisher  Information  matrix, 
see  Ref.  3).  If  we  define 

B(t)  =  3(y/m<y/£)  (2.10) 


as  the  information  matrix  for  a  single  observation  at  time  t,  then  Eqs.  (2.  6)  or  (2.  5) 


*  It  is  also  called  the  normal  matrix. 
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shows  that,  for  independent  observation  errors,  the  information  matrices  add.  Thus 
if  we  make  two  sets  of  observations  with  information  matrices  Bj  and  B ^  and  if  the  two 
observation  noise  processes  are  mutually  independent,  then  the  total  htformation  matrix 
B,  is  simply 


B  =  B 


1 


(2.11) 


Often  knowledge  of  the  value  of  £0  is  obtained  from  sources  other  than  just  the 
observations.  If  we  have  such  aprtcri  information,  it  can  often  be  modeled  as  a  set 
of  values 


_  .  o 

—cl  priori  —  -ki pHori 


(2. 12) 


where  r?  .  .  is  a  zero  mean,  Gaussian  random  vector  with  covariance  matrix 

- opr  tort 

B  *  .  .  If  »?  ..  .  i  is  independent  of  the  observation  noise,  the  additive  property 

cl  priori  —cl  prior  t 

of  information  makes  B  +  B  .the  information  matrix  corresponding  to  both  the 

cl  p  non 

a  priori  information  and  the  observations. 

Our  results  are  not  completely  tied  to  a  Gaussian  assumption.  Using  the 
Taylor  expansion  approach  of  Eq.  (2. 7),  the  covariance  matrix  of  Eq.  (2.  9)  is  that 
obtained  by  least  squares  unbiased  estimators  (minimum  variance  linear  estimators) 
for  any  "white"  but  non-Gaussian  r?(t).  However,  for  simplicity,  we  continue  to  employ 
the  Gaussian  assumptions. 

Henceforth  we  call  the  matrix  2  defined  by  Eq.  (2.  9)  and  Eq.  (2. 6)  the  covariance 
matrix  of  the  estimate  errors.  Our  preceding  discussions  show  this  to  be  a  valid 
statement  provided  1)  we  use  maximum  likelihood  (matched  filter)  estimation  techniques 
and  2)  the  signal -to-noise  ratio  and  amount  of  aprtoH  information  are  large  enough  to 
give  small  errors  in  the  estimates  of  the  parameters.  If  either  condition  is  violated, 
our  results  are  actually  a  bound  on  the  obtainable  accuracy. 
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3.  SINGLE  TARGET  INFORMATION  MATRIX 

The  formulae  of  Sec.  2  are  now  applied  to  a  general  model  for  a  single  target. 

A  radar  transmits  a  modulated  sine  wave  with  carrier  frequency  and  time  duration 
T.  Let  r(t)  denote  the  range  between  the  radar  and  a  moving  point  reflector  at  time  t. 
Let  6(t)  denote  the  time  delay  between  transmission  and  receipt  of  the  signal;  to  be 
explicit,  a  signal  received  at  time  t,  was  transmitted  at  time  t-6(t).  r(t)  and  6(t)  are 
related  by 

c«(t)  =  2r(t (3.1) 


where  c  is  the  speed  of  light.  Assume  the  equations  of  motion  of  the  reflector  are 
known  to  within  the  values  of  the  unknown  parameters.  By  •  •  •  ,  0^  .  Let  B_  be  the 
column  vector 


e 


i 


0. 


(3.2) 


The  time  delay  at  time  t  is  written  as  6(t,0).  We  do  not  specialize  our  results  to  any 
specific  form  for  6(t,0).  However,  a  common  model  for  accelerating  targets  is 


<5(t,0)  =e1  +02(t-T)  + 


(3.3) 


where  t  is  some  chosen  reference  time.  For  small  accelerations,  0^,  6^  and  0^  are 
approximately  proportional  to  the  target's  range,  range  rate  and  range  acceleration 
at  time  r. 

Let  oi  denote  the  amplitude  and  j3  the  phase  of  the  reflected  signal  as  received 
at  the  transmitter.  Both  ot  and  (3  are  considered  unknown  parameters.  Let  £  denote 
the  p  +  2  dimensional  column  vector  containing  the  elements  of  0  and  a  and  /3.  Let£° 
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and  thus  0°,  cP  and  /3°  denote  the  true  but  unknown  values  of  these  parameters. 

The  carrier  is  modulated  by  both  amplitude  and  frequency  modulation.  Let 
v(t)  denote  the  amplitude  modulation  and  let  u(t)  denote  the  frequency  modulation. 

Assume  both  quadrature  components,  z^(t)  k  =  1,  2,  of  the  received  signal  are 
observed  in  the  presence  of  additive,  broad  band,  zero  mean  Gaussian  noise.  Then* 

zk(t)  =  yk(t,f0)  +  nk(t)  k  =  1,2  0SI5T  (3.4) 

where  rj  (t)  and  rj  (t)  are  zero  mean,  stationary,  white,  mutually  independent,  Gaussian 

1  1  2 
stochastic  processes  with  two-sided  noise  power,  a  and 

yk(t,f)  =  orv(t,£)  cos[o>M)+  fi  -  ]  k  =  1,2.  (3.5) 

The  exact  form  of  v(t,£)  and  w(t,£)  depends  on  the  explicit  problem.  If  both  phase  and 
envelope  characteristics  of  the  observed  signal  are  used  to  form  the  desired  estimates, 

v(t,0)  =  v[t-6(t,0)]  (3.6) 

pt-6(t,0) 

w(t,0)  =  \  +  u(T)]dT  .  (3.7) 

0 

In  Sec.  5,  the  above  and  forms  corresponding  to  other  cases  will  be  discussed. 

Equations  (3.4)  and  (3.  5)  define  a  nonlinear  parameter  estimation  model  of  the 

type  discussed  in  Sec.  2.  Therefore,  2  ,  the  covariance  matrix  of  the  errors  made  in 

o  « 

estimating  from  observing  z^(t)  0  ^  t  ^  T  is  given  in  the  sense  of  Sec.  2  by 


*  The  time  interval  for  which  y^(t,  j[°)  is  nonzero  is  really  t^(£°)  ^  t  ^  T  +  6(T,0°) 

where  tj(0°)  is  the  first  time  for  which  t-6(t,£°)  is  nonnegative.  However,  for  most 
problems  of  interest,  an  interval  0  ^t  <  T  is  a  satisfactory  model  for  analysis.  If  not 
the  necessary  modifications  are  obvious. 
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(3.8) 


Yt)=Bl,£(t)  +  B2^> 


(3.9) 


(3.10) 


BkYt)=  — 


k  =  1,  2 


(3.11) 


where  3(yk/£)  denotes  the  column  vector  of  partial  derivatives  of  yk(t,  f)  wit^ 


respect  to  the  components  of  £,  evaluated  at  £  =  £°.  B  (t)  is  the  information  obtained 

o 

at  time  t  from  observing  both  quadrature  components  of  the  received  signal  while 

is  the  total  information  obtained  over  the  time  interval  O^t^T.  B„  .  . 

£,  a  prt  ort 

is  the  a  priori  information.  Our  final  interest  is  in  2^  but  we  often  confine 
discussions  to  just  B^(t)  as  it  is  the  basic  quantity. 

If  we  partition  £  as 


then 


(3.12) 


3(yk/£)  = 


?f(yk/a) 

s<yk/0> 


3(yk/£) 


(3.  13) 


where  B^y^/cv)  and  %(y^/P)  are  scalars  while  I5(y  /0_)  is  a  p-dimensional  column 
vector.  Here  and  in  the  rest  of  the  report,  we  use  dotted  lines  to  indicate  the  partitioning 
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of  vectors  and  matrices  (which  is  not  always  as  done  in  Eqs.(3. 12)  and  (3. 13).  Define 


nk  =  u>(t,£)  +  /3-(k-l)7r/2.  (3.14) 

Then  from  Eq.  (3.  5) 

®(yk/a)  =  v(t,0)  cos 
5(yk/0)  =  “tw(t,£)  sin 

S(yk/£)  =  «[3(v/£)  cos  n  -  v(t,£)t!(o;/£)sin  B  ]  . 


Combining  Eqs.  (3.  10)  and  (3.  11)  with  Eqs.  (3. 15),  (3. 16)  and  (3.  17)  gives  after 
manipulation  and  cancellation 


v2(t) 

2 

a 

0 

1 

1  ~v(t)t5'(v/£) 

1 

v2(t) 

v2(t)3'(u>/£) 

2 

(t)  =  — 

—  a 

—  —  - 

-  —  —  —  —  — 

1 

! 

1 

1 

j  t5(v/£)  tT'(v/£)  +  v2(t)  ?(w/£)  3#(w/£) 

(3.15) 

(3. 16) 
(3.  17) 


where  Eq.(3.  18)  is  partitioned  corresponding  to  Eq.(3.  12)  with  the  lower  right  hand 
element  a  p  x  p  matrix.  We  indicate  only  half  of  the  off-diagonal  terms  as  V0  is 
a  symmetric  matrix.  Equation  (3.  18)  is  to  be  evaluated  for  £  =  £  (i.  e. , 
a  =  a  ,  £  =£°)  but  for  simplicity  we  do  not  include  this  or  the  dependence  of  v(t)  on  £ 
in  our  notation. 

In  the  following  sections  we  work  solely  with  the  time  delay  6(t,0)  and  in 
particular,  ?(6/£).  If  the  corresponding  range,  r(t,£)  is  given  instead,  Eq.(3.1)must 
be  used.  A  useful  formulae  obtained  from  Eq.(3.  1)  is 
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[c  +  r(t  -5(t, £_),£)  ]S(6/0)  =  2&(r/0) 


where 


r(t  -6(t,0),0) 
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4.  EFFECTS  OF  UNKNOWN  PHASE  AND  AMPLITUDE  ON  TARGET  PARAMETERS 


Our  analysis  includes  the  variances  of  estimates  of  the  phase,  0  and  the 
amplitude  a.  These  are  sometimes  of  interest;  for  example,  the  variance  of  an 
estimate  of  a  is  a  measure  of  detection  ability.  However,  for  many  accuracy  studies 
a  and  /3  are  nuisance  parameters  and  we  are  really  interested  in  2=13*,  the 

V  U 

covariance  matrix  of  the  errors  made  in  estimating  the  target  parameters  0.  B  can 

♦  _ 

be  obtained  by  partitioning  as 


and 


I 


B11  I 

B12 

1 

1 

_  1 

fill  ' 
2  1 

_  | 

s12 

1 

__  1 

f  CN 

1  “ 

1 

(4.1) 


(4.2) 


If  the  partitioning  is  done  as  indicated  by  the  dotted  lines  in  Eq.(3. 18), 


2 


22 


(4.3) 


The  formulae  for  inverting  partitioned  matrices  give 


a:2V1-B22-B2y,)-,B12 


(4.4) 


Applying  Eq.  (4.  4)  to  Eq.  (3. 18)  gives 


*  If  a  vrtort  information  is  present,  we  partition  B,_  +  B..  .  . 

£  ^apnori 


11 


(4.5) 


where 


—  B0=  f  [g(v/e)S/(v/£)  +  v2(t)3(w/£)S/(w/£)  dt 

rv  —  ^  L  J 


1 

_  (-*!  .  /-<  (1  / 

f  v2(t)dt 

U  U  T  u  VJ 

L  V  V  co  co  _ 

0 


G 

v 


I 


T 

v(t)  t5(v/£)dt 


0 


=  y  v2(t)O(w/0)  dt 
0 


(4.6) 


(4.7) 


The  matrix,  G  G  '  +  G  G'  is  nonnegative  definite  and  measures  the  loss  in 
V  V  U)  co 

information  about  6_  that  results  from  an  unknown  phase  and  amplitude.  As  we  shall 
see,  this  information  loss  is  zero  in  some  special  cases. 

We  do  not  use  B  in  our  further  development  as  consideration  of  the  full  B 

_  L 

matrix  leads  to  much  simpler  results  and  a  better  understanding  of  the  problem. 
However,  Eq.(4.  4)  can  be  employed  at  any  point  B^  is  of  interest. 


12 


5. 


SINGLE  TARGET  EQUIVALENT  LINEAR  MODELS 


We  now  discuss  linear  equivalent  models  for  three  cases  of  a  single  target: 

R-F  Phase  Used  for  all  Parameters 
R-F  Phase  Used  for  Motion  Parameters 
R-F  Phase  Not  Used. 

These  cases  do  not  exhaust  the  range  of  possibilities  but  are  sufficiently  illustrative 
to  indicate  how  other  situations  are  handled.  The  choice  of  which  case  to  use  naturally 
depends  on  the  problem  under  investigation;  in  particular,  the  type  of  the  equipment 
available  and  the  nature  of  the  ambiguity  problem. 

In  Sec.  2,  we  discussed  the  linearization  of  a  nonlinear  problem  by  truncation 
of  a  Taylor  series  expansion.  In  this  sense,  we  could  form  an  equivalent  linear  model 
for  our  radar  problem  by  using  the  partial  derivatives  of  Eqs.  (3. 15),  (3.  16)  and  (3. 17). 
However,  we  follow  a  more  productive  course  and  hypothesize  two  sets  of  observations 


za(t)  =  2a'(t)I°  +  VC> 

(5.1) 

zf(t)  =  2[(t)£°  +  T)f  (t) 

(5.2) 

where  n  (t)  and  t]  ( t)  are  zero  mean,  stationary,  white,  mutually  independent  Gaussian 
a  i  ,, 

2  o 

stochastic  processes  with  two-sided  noise  power  o  ,  £  is  the  parameter  vector  to  be 
estimated,  and  £a(t)  and  £^.(t)  are  column  vector  functions  of  time.  We  know  from 
Sec.  2  (or  Eqs.  (3.  8)  through  (3.  11))  that  the  information  obtained  from  observing  z  (t) 
and  7^(t)  at  time  t  is, 


Y°  *  Ba,  C<'>  +  Bf, 


(5.3) 


where 


Bf  £  (t)  =  -y2f(t)2'(t)  . 
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Equations  (5. 1)  and  (5. 2)  are  an  equivalent  linear  model  for  a  radar  problem  when  the 
vector  functions  ^(t)  and  ^(t)  are  chosen  to  make  B^(t)  of  Eq.  (5. 3)  the  same  as  the 
B^(t)  of  the  original  radar  problem. 

In  a  certain  sense,  we  choose  £^(t)  and  ^(t)  such  that  Ba  ^(t)  represents  the 
information  obtained  from  the  amplitude  modulation  while  B  (t)  represents  the 

M  b 

information  obtained  from  the  frequency  modulation  and  the  R-F  phase.  As  will  be 
seen,  this  dichotomy  is  not  precise  as  ^(t)  is  dependent  on  the  amplitude  modulation 
v(t)  as  well  as  the  frequency  modulation,  u(t).  However  for  many  situations,  v(t)  itself 
determines  the  energy  content  of  the  signal  while  it  is  the  structure  of  ^-v(t)  that 
determines  the  effectiveness  of  the  amplitude  modulation.  Admittedly,  this  point  of  view 
is  not  universally  valid  but  in  the  author's  opinion,  it  is  of  sufficient  value  to  call  ^(t) 
and  gf(t)  equivalent  models  for  amplitude  and  frequency  respectively. 

If  any  a  priori  information  exists,  the  equivalent  linear  model  of  Eqs.  (5. 1) 
and  (5.2)  is  supplemented  by  the  additional  observations  (see  Eq.  (2. 12) 

_  t°  , 

-a  pri  ori  -  -a priori  (5.4) 


with  B„  .  .  the  corresponding  information.  For  example,  if  . 

f,a  priori  y  p  0,c  priori 

denotes  the  covariance  matrix  of  an  estimate  of  £  made  from  previous  tracking  of  the 
target  and  if  we  partition  as 


t  = 


(5.5) 


then' 


.-1 


*  Note  that  B  ‘  J  does  not  exist.  In  the  same  sense,  T.n  .  M  does  not  have 
t,  a  priori  6_,  a  priori 

to  exist. 
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0 


B 


La  priori 


0 

0 

•  o 

1 

1 

1 

1 

0 


-1 


'0,  a  priori 


Previous  tracking  might  also  provide  useful  information  on  the  target  cross  section 

which  implies  a  priori  information  on  the  amplitude  at.  A  priori ■  information  can  also 

arise  from  other  sources  such  as  equipment  which  maintains  a  measure  of  the  phase. 

2 

If  a.  .  .  denotes  the  variance  of  these  equipment  errors  and  other  phase 

p  t  CL  Jp7*u07*u 

uncertainties,  then 


B 


a  priori 


0 


0 


I 


0,  a  priori 


0 


0 


0 


(5,6) 


is  the  corresponding  apriori  information. 

The  limiting  case  where  the  carrier  frequency,  w  ,  is  much  larger  than  the 
effective  signal  bandwidth  is  of  special  interest.  This  can  be  obtained  by  simply 
dropping  the  lower  order  terms  in  that  affect  the  information  on  0_.  It  is  difficult 
to  specify  the  range  of  carrier  frequency  to  bandwidth  for  which  this  approximation  is 
useful  as,  in  general,  the  nature  of  the  target's  motion  is  also  a  contributing  factor. 
However,  the  large  c approximation  is  very  useful  for  a  wide  range  of  problems. 
The  large  case  can  also  be  obtained  by  appropriately  modifying  the  definitions 
of  v(t,£)  and  ta(t,0_). 
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# 


To  illustrate  the  equivalent  linear  model,  we  begin  with  the  general  case,  It 
is  seen  by  inspection,  that  for 


*a(,)  = 


£f(t)  =  <*v(t) 


the  B  (t)  obtained  from  Eq.(5.3)  is  equal  to  the  B  (t)  of  Eq.  (3.  18), 

e>  b 

We  now  consider  three  special  cases  and  obtain  explicit  formulae  for 
3(w/£)  and  $5(v/0)  . 
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R-F  Phase  Used  For  All  Parameters 


For  the  case  when  the  R-F  phase  is  used  for  all  parameters, the  equations 
for  v(t,0)  and  ci>(t,0)  are  given  by  Eq.(3.6)  and  (3.  7).  The  resulting  partial  derivatives 
are 


S(v/£)  =  -v(t)  t5<6/0)  (5.7) 

3(aV0)  =  -[a>c  +  u(t)]  3(6/0)  (S.  8) 


where 


v(t)  = 


dv(t) 

dt 


The  time  functions  of  Eqs.  (5.  7)  and  (5.  8)  are  actually  evaluated  at  time  t-6(t,  0°)  but 
evaluation  at  time  t  is  satisfactory  for  most  cases.  Therefore  we  simplify  our  notation 
by  dropping  the  0°  dependence.  Evaluation  of  Eq.  (3.  18)  for  Eqs.  (5.  7)  and  (5.  8)  gives 


a 

V»*T 

-  a 


v2(t) 

2 

a 


0 


v2(t) 


-V(t)t5'(6/0) 

cv 


v  (t)Cwc  +  u(t)]3#(6/0) 


|v2(t)  +v2(t)[a.c+u(t)]^  W/m'( 6/9) 


(5.9) 


where  the  partitioning  is  done  as  in  Eq.(5.  5).  By  inspection  we  see  that 
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* 


v(t) 

Of 

0 

v(t)  3(6/0) 

fn 

<£.f(t)  =  ttv(t) 


1 


(5.  10) 


(5.11) 


[wc  +  u(t)]  3(6/0) 


combine  with  Eq.(5.3)  to  give  the  same  B^(t)  as  that  of  Eq.  (5.  9).  When  the  R-F  phase 
is  used  for  all  parameters,  a  priori  information  on  the  phase  is  essential  and  this  is 
m<~.  leled  as  in  Eq.  (5.6).  Thus  our  equivalent  linear  system  consists  of  Eqs.  (5.  1)  and 
(5.2)  with  the  <£_JX)  and  g?  (t)  of  Eqs.  (5.  10)  and  (5.  11),  combined  with  the  a  priori 
model  of  Eqs.  (5.  4)  and  (5.  6). 

2 

If  the  phase  ft  is  known  perfectly,  i.e. ,  rr  .  .=  0,  the  row  and  column 

7  ft,  a  priori 

of  Eq.  (5.  9)  corresponding  to  ft  are  simply  c  Meted.  If  in  addition  we  use  the  large 
approximation,  our  equivalent  linear  model  for  estimating  0_  is  simply 


where 


zf(t)  =  £f'(0£°+  tjf(t) 
<2_f(t)=  av(t)w  3(6/0) 


(5. 12) 


t  =e 
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R-F  Phase  Used  For  Motion  Parameters 


Consider  the  case 


6(1,0)  =0+6  (1,0  ) 

-  1  ni  -m 


(5.  13) 


where 


»,  =0<t.0)/t=T 


(5.  14) 


and  0  is  the  p-1  dimensional  column  vector  containing  all  elements  of  0  except  0  , 

— m  — 

is  the  time  delay  at  time  t  and  6  (t,£  )  describes  the  change  in  6(t,0_)  with  the 
motion  of  the  target.  Thus  the  elements  of  0_^  are  the  parameters  of  the  target's 
motion.  For  the  accelerating  target  example  of  Eq.(3.3), 


6  (t,0  )  =  0  (t-T)  + 

m  —  m  2. 


0  3(t-r)2 
2 


0 


— m 


In  many  situations,  lack  of  knowledge  of  phase  /3  or  more  often  the  ambiguity  problem 

prevent  the  effective  use  of  the  R-F  phase  in  the  determination  of  0^  although  it  can  be 

used  for  the  target  motion  parameters  0  .  For  such  cases  we  write 

— m 


v(t,  0)  =  v(t-0  -6  (t,  0  )) 

—  l  m  — m 


(5.  15) 


uo(t,  0)  =  00  t  -00  <5  (t,  0  )  + 
—  c  cm  — m 


^t  -0-6  (c  ,0  ) 

f  1  m  "mu(s)ds. 


(5.  16) 


0 


In  Eq.  (5.  lb)  the  oo^O  ^  term  has  been  combined  with  the  phase  shift  (1  (see  Kq.  (3.  5))  ; 
that  is,  +  w^0  ^  1S  considered  the  unknown  phase  shift  rather  that  just  /5 .  This 
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change  of  variables  is  made  for  convenience,  not  necessity.  For  example,  if 
Eq.  (5.  13)  is  substituted  into  Eq.  (5.  9)  and  the  partitioning  technique  of  Sec.  4  applied 


to  the  resulting  B^,  the  same  2^  we  obtain  by  the  change  of  variables  results  as  many 
terms  cancel.  However,  if  we  keep  only  the  highest  order  terms  in  the  B^ 
obtained  from  Eq.  (5.  9)  becomes  singular  for  fift.OJ)  of  Eq.  (5.  13).  The  change  of 
variables  side  tracks  these  unnecessary  complications. 

If  we  partition  0  as 


0 

-m 


then  from  Eqs.  (5.  15)  and  (5.  16) 


?[(v/fl)=  -v(0 


m  /0  ) 

m  —  m 


(5.17) 


u(t) 


ff(w/£)  =  - 


o  +u(t)]?(6  /O  ) 
c  m  —  m 


(5.18) 


The  resulting  B^(t)  of  Eq.  (3.  18)  is* 


*  A  constant  frequency  modulation  u(t)  =  k,  is  equivalent  to  a  carrier  frequency  of 
u;  +  k  with  u(t)  =  0.  This  obvious  physical  deduction  is  not  evident  from  Eq.  (5.  19) 

but  use  of  the  matrix  partitioning  of  Sec.  4  shows  it  to  be  true  as  various  terms  cancel. 
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(5.  19) 


.  01 

vt>=- 

~  a 


V(t>  0 

cy 


-£v(t)v<t) 


v2(t)  -v2(t)u(t) 


.  2.  .  2,  .  2.  v 

v  (t)+v  (t)u  (t) 


J 


-v(t)v(t)3'(6  /0  ) 
a  m  -m 


-{v2(t)[cu  +u(t)]}ff/(fim/0  ) 
c  m  - m 


{v2(t)+v2(t)u(t)[wc+u(t)]}?:'(6m/£m) 


{v2(t)+v2(t)[w  +u(t)]2}  W/en'(6  /o) 
c  m  —  m  m  -m 


where  we  have  partitioned  as 


t  = 


a 

P 
0 . 


0 

-m 


By  inspection  we  see  that  the  corresponding  forms  of  Eqs.  (5.  1)  and  (5.  2)  are 


0  (t) 
—a 


=  O' 


-v(t) 

O' 


0 

v(t) 


v(t)  m/e)  ( 

I _  m  — m_| 
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<0(t)  =  «v(t) 

— I 


0 


u(t) 


L 


O  +  u(t)]3(6/(!m) 

c  m  —  m 


For  the  large  w 


c 


approximation, 


-v(t) 

Of 

0 

v(t) 

0 


=  «v(t) 


0 

1 

u(t) 


W  ?(6  /0  ) 

c  m  — m 


As  mentioned  earlier.the  large  w  case  can  also  be  derived  directly  by  changing  the 
definitions  of  v(t,£)  and  o;(t,£).  For  example,  the  present  case  can  also  be  obtained 
using 
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v(t ,o)  =  v (t  -n 

ft_0i 

u(t.O)  =  ca  t  -  co  6  (t  ,0  )  +  \  u(s)ds 

-  c  cm  — m  J 

0 


and  substituting  the  resulting  n(v/0_)  and  %(<jo/0)  directly  into  Hq.(3.  18).  Reference  2 
uses  this  approach. 

R-F  Phase  Not  Used 

If  the  R-F  phase  is  not  used  at  all,  the  phase  parameter  ft  is  dropped  from 
our  model  by  removing  the  corresponding  row  and  column  of  Eq.  (3.  18).  Using 


v(t ,0)  -■  v[t-6(t,  £)] 

pt-^(t.O) 

w(t,  0)  =  w  t  +  \  u(r)dT 

c  J 
0 


and  partitioning  as 


we  obtain 


— 

1 

v2(t) 

2 

a 

1 

1 

i 

i 

-^v(t)v(O?5/(6/0) 

l 

! 

1 

l 

[v2<0  +  v2(t)u2(t)]3(6/e)  !5'(6/i> 

(5.20) 
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The  corresponding  linear  model  is 


v(t) 

O' 


a 


v(t)?(6/£) 


<£f(t)  =  av(t) 


u(t)r!(6/0) 
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6.  MULTIPLE  TARGETS 


In  many  applications,  target  resolution  rather  than  single  target  accuracy  is  of 
prime  interest.  Target  resolution  refers  to  the  ability  of  the  radar  to  deter  mine 
(resolve)  the  properties  of  separate  targets  when  it  observes  the  sum  of  the  reflections 
from  each  target.  One  approach  to  this  problem  is  the  inspection  of  the  ambiguity 
function.  However,  a  more  precise  technique  is  to  simply  extend  the  single  target 
analysis  to  the  multiple  target  case  and  thereby  investigate  the  following  question: 

Given  a  noise  corrupted  observation  of  the  sum  of  the  signals  reflected  from  several 
targets,  how  well  can  the  parameters  of  all  the  targets  be  estimated?*  (Reference  7 
uses  this  technique  to  investigate  angular  resolution.) 

For  simplicity,  we  restrict  our  explicit  equations  to  the  case  of  just  two  targets. 
The  general  multiple  target  case  is  conceptually  a  trivial  extension  but  the  notational 
complications  tend  to  obscure  the  basic  ideas.  In  the  two-target  case,  there  are 
2(p+2)  unknown  parameters  which  we  denote  by  the  column  vector 


and  of  course,  a  ,  f}  and  0  denote  the  amplitude,  phase  and  range  parameters  of  the 

th  J  J  J 

j  target. 

*  Our  single  target  results  can  be  interpreted  as  the  curvature  of  the  ambiguity  ianction 
at  the  true  values  of  the  parameters  (see  Ref.  2).  Thus  the  multiple  target  analysis  can 
be  considered  the  curvature  of  a  higher  dimensional  ambiguity  function. 
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For  the  two-target  case, the  two  quadrature  components  corresponding  to  Eq.  (3.4) 
are  given  by 


Zk(0  =  Yk(t,Z  )  +  nk(t)  k  =  1,2 


where 


Yk(uz)  =  yk(ui,)  +  yk(t-l2) 


and  yk(t,  f)  is  given  by  Eq.  (3.  5).  It  follows  that,  partitioning  as  in  Eq.  (6.  1), 


?(Yk/Z)  = 


3<yk/I,) 


Wyk/I2> 


(6.3) 


where  3(y  /£ .)  is  the  column  vector  of  partial  derivatives  of  y  (t,£.)  with  respect 
k  j  k  j 

the  elements  of  Using  Eq.  (3.  11)  and  (3.  10)  we  have 


to 


B 7(t)  =  — j-  {?(Y1/Z)3'(Y1/Z)+3(Y2/Z)3'(Y2/Z)}  . 


Substituting  Eq.  (6.  3)  gives 


BM 


(t) 

1  C(t) 

1 

1 

1 

t,; 

—2 

i 

(6.  4) 


where  partitioning  is  as  in  Eq.  (6.  1),  B  (t)  and  B  (t)  are  the  information  matrices  which 

-1  —2 

would  arise  if  the  targets  were  observed  separately,  and 
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r 


9 


C(t)  = 


(6.5) 


Using  Eq.  (3.  13), 

8(V'V'5(yk/0!2>  ’5<yk/'V8(yk/,,2) 

I 

3(yk//3,)8(yk/»2)  3<yk/(5,)3(yk//s2)  I 


3(yk/a2)3(yk/e ,)  3<yk//3  2)%k/i ,) 


3(yk/“|)3'(yk/£2) 

3(yk//s1)3'(yk/£2) 


(6.  6) 


3(yk/£,)3'(yk/92) 


where  partitioning  is  as  in  Eq.  (6.  2).  The  necessary  partial  derivatives  are  given  by 
Eq.  (3.  15)  through  Eq.  (3.  17)  with  the  obvious  subscript  additions.  Define 

A(t)  =  w(t,0  {)  -  w(t ,02)  +  (3[  -  p2  •  (6.7) 


For  simplicity  we  do  not  indicate  the  dependence  of  A(t)  on  p^  and  p^ 

define 


Also 


v  (t)  =  v (t .  0 . )  j  =  1,2  . 

J  “J 

Substituting  Eqs.  (3.  15),  (i.  16),  (3.  17)  into  Eq.  (6.  6)  then  gives  after  reduction  by 
use  of  trigonometric  identities, 
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28 


Eq.  (6.  8) 


Thus  the  information  matrix,  B  (t)  of  Eq.  (6.  4)  is  completely  specified  by  lup .  ( t> .  8)  and 

L/ 

(3.  18). 

As  in  the  single  target  case,  we  have 


(B„  +  B„  .  .  ) 

Z_  Z,  a  priori 


-1 


fTt 

\  B  (t)cit 
0 


where  T7  is  the  covariance  matrix  of  the  errors  in  the  unbiased  estimate  of  Z  .  In  the 
multiple  target  case,  the  possible  range  differences  spread  out  the  time  duration  of 
observations  and  this  is  indicated  by  making  the  observation  interval  0  to  T  where 
is  simply  chosen  large  enough  to  encompass  the  actual  observation  interval. 

One  important  property  can  be  observed  directly  from  Eq.  (6.  4).  Define,  with 
partitioning  as  in  Eq.  (6.  J) 


A  = 

L 

where  with  partitioning  as  in  Eq.  (6.  2) 


i 


i 


A 


2 


A 

J 


i 

1  0 

1 

1 

O' 

1 

j 

1 

_l_ 

l 

|  O'  1 

j 

_ 

1 

where  1  is  a  p  by  p  unit  matrix.  Eq.(6.4)  can  be  written  as 


< 
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Bz(t)  =  A  Bz(t)  A 

A 

where  B„(t)  does  not  depend  on  the  amplitude  parameters  and  Ct  Thus  if  there 
is  no  aprto^t  information, 


which  shows  that  the  amplitude  (cross  section),  0t^v  of  the  second  target  has  no  effect 
on  the  accuracy  with  which  the  parameters  of  the  first  target  can  be  estimated  and  vice 
versa.  *  (See  also  the  angular  resolution  problem  of  Ref.  7.)  This  property  is  a 
result  of  our  assumption  that  observation  noise  is  the  only  source  of  error.  In  practice, 
uncertainties  in  waveform  shape  and  imperfections  in  tht  data  processing  equipment 
will  introduce  other  errors  which,  for  large  cross  section  targets,  may  not  be  negligible. 

Now  let  us  consider  an  equivalent  linear  model  for  the  two-target  case.  It  is 
convenient  to  resort  to  a  vector  formulation.  Assume  the  q^ -dimensional  column  vectoi 
£(t)  is  observed  where 

z(t)  =  *  \t)Z°  +  Mt)  0£tsTt  (6. 9) 

where  Z°  is  a  q^-dimensional  column  vector  of  unknown  parameters,  4>(t)  is  a 

^1  by  ^2  matr*x  known  time  functions  and  rj( t)  is  a  q(  dimensional,  zero  mean  white 

Gaussian  noise  process  with 

E[n(t)9'(t-T)]=  <5(r)  I  (6.10) 

where  in  Eq.  (6.  10),  6(t)  is  the  Dirac  delta  function  and  not  a  time  delay.  For  the 
model  of  Eq.(6.  9)  and  (6.  10)  we  have,  (for  the  corresponding  discrete  time  case,  see 
Ref.  6) 

*  Unless,  of  course,  0^  is  known  to  be  zero. 

.10 


(0.11) 


(0.12) 


(6.13) 


The  equivalent  linear  model  for  the  single  target,  Eqs.  (5.  1)  and  (5.  2),  can  be  written 
in  terms  of  this  vector  model  as 


(6.14) 


(6.15) 


(6.  10) 


(6.17) 
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For  the  two-target  case,  we  define,  b>  analogy  with  the  s  ngle  target  asc  of 


Sec.  5. 


v.(t) 

J _ 

a. 

J 


S(v/0  ) 


(6.  18) 


(6.19) 


We  further  define 


e'(t)  = 


£ 


f.  1  <0 


r(t)  = 


0 


1 

0 

cosA(t) 

sinA(t) 


0 


0  <,.2(t) 


£f,  2(t) 


0 

1 

-sinA(t) 

cosA(t) 


(6.20) 


(6.21) 
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Then  by  algebraic  manipulation,  it  can  be  shown  that  B  (t)  of  Eq.  (6.  4)  can  be  wr  itten 

Lj 


B  7(t) 


~ 7~  0(t)H(t)  1 
a 


0 


'(t)6'(t)  . 


1 


Comparison  of  Eq.  (6.  22)  with  Eq.  (6.  13)  shows  that  two  different  linear 
models  are  of  interest.  For  both  models  we  have 


q2  =  2  (p  +  2) 


One  possible  linear  equivalent  model  is  then  given  by 


and 

$'(t)  =  H'(t)9'(t) 


or 


/ 

£ 


a,  1 


(t) 


*'(t)  = 


cosA(t)  (£_’  (t)  +  sinA(t)  q>'  (t) 


-sinA(t)  ^  2(t)  +  cosA(t)  2(t) 


The  second  linear  equivalent  model  is  given  by 


(6.  22) 


equivalent 


(6.  23) 


(6.24) 
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ql  =  4 


*'<t)  =  o'(t) 


(6  25) 


and 


or 


<  -  -4 

—  a 


1 

0 

cosA(t) 

sinA(t) 


“I  1  V 

v  =  — 

a 


0 

1 

-sinA(t) 

cosA(t) 
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(6.26) 


The  equivalent  linear  model  of  Eq.  (6.  23)  and  (6.  24)  differs  from  the  single  target 
case  due  to  the  extra  parameters  and  a  more  complicated  expression  for  <J>'(t).  The 
equivalent  linear  model  of  Eq.  (6.  25)  and  (6.  26)  employs  four  instead  of  two  observations 

and  has  the  following  interesting  property:  If  of  Eq.  (6.  26)  were  — y-  I,  the  model 

—  a 

would  correspond  to  the  case  wherein  the  reflections  of  the  two  targets  are  observed 
separately.  Thus  all  the  effects  of  the  multiple  target  are  contained  in  the  cross 
correlation  terms  of  Eq.(6.26).  A  disadvantage  of  the  second  model  is  that  E’1  of 
Eq.  (6.26)  is  a  singular  matrix;  that  is,  does  not  exist. 

Since  the  multiple  target  model  employs  the  same  <p  (t)  and  q> A t)  used  for  the 
single  target  case,  the  discussions  in  Sec.  5  on  the  various  special  cases  need  not  be 
repeated.  The  additional  term,  A(t),  of  Eq.  (6.7)  is  easily  found  from  the  equations 
ce(t,£)  given  in  Sec.  5. 

In  a  multiple  target  environment,  the  question  of  how  many  targets  are  present 
is  very  important.  For  exa mple, if  t\so  targets  are  assumed  when  actually  there  is  only 
one,  the  estimates  of  the  single  target's  parameters  are  degraded.  The  variances  of 
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the  estimates  of  the  amplitude  parameters,  »  ,  measure  how  well  the  number  of  targets 
actually  present  can  be  decided.  (The  corresponding  discussion  in  Kef.  7  is  more 
detailed. ) 

The  matrix  partitioning  t  scussed  in  Sec.  4  is  both  applicable  and 

useful  in  the  multiple  target  case.  For  example,  all  the  parameters  of  a  second  target 
might  be  considered  nuisance  parameters. 

In  many  single  target  problems,  the  actual  values  of  the  unknown  parameters 
do  not  affect  the  information  matrix  However,  for  multiple  targets  the  information 
matrix  and  hence  the  covariance  matrix  is  strongly  dependent  on  the  relative  values  of 
the  parameters  of  the  different  targets.  Thus  in  order  to  obtain  a  true  picture  of  target 
resolution  ability,  it  is  necessary  to  investigate  the  effects  of  changes  in  target  separation. 
(Reference  7  indicates  the  type  of  results  which  may  be  expected.) 


7.  DISCUSSION 


We  have  used  the  information  matrix  as  a  basis  for  investigating  the  accuracy 
and  resolution  ability  of  a  wide  class  of  radar  signals.  The  inverse  of  this  information 
matrix  is  a  lower  bound  on  the  covariance  matrix  of  the  errors  obtainable  from  any 
unbiased  estimate  of  the  unknown  parameters.  In  theory,  this  bound  is  approached  as 
the  signal-to-noise  ratio  increases  if  maximum  likelihood  parameter  estimation 
procedures  are  used.  However,  for  very  large  signal-to-noise  ratios,  errors  in 
factors  such  as  waveform  shape  and  data  processing  equipment  become  as  important 
as  those  due  to  observation  noise.  Thus  this  bound  is  obtainable  only  for  a  limited 
range  of  signal-to-noise  ratios.  We  have  not  attempted  to  discuss  methods  of  imple¬ 
menting  the  required  estimation  procedures  as  there  are  many  possibilities  and  the 
choice  of  technique  depends  very  much  on  the  explicit  problem  being  investigated. 

The  inverse  of  the  information  matrix  is  usually  of  prime  interest.  In  some 
special  cases  an  analytic  matrix  inversion  is  practical  (see  for  example  Refs.  2  and  8) 
but  there  seems  to  be  little  hope  for  the  general  models  we  have  considered.  Thus  if 
explicit  closed  form  expressions  are  desired,  the  problem  of  interest  must  be 
restricted  further  and  a  potentially  very  difficult  analysis  attempted.  The  matrix 
partitioning  technique  of  Sec.  4  may  be  helpful. 

In  the  single  target  case,  the  equivalent  linear  model  is  a  convenient  starting 
point  for  the  design  of  modulations  which  are  optimized  with  respect  to  estimate 
variances.  *  A  general  design  theory  for  such  problems  is  discussed  in  Ref.  9  and  some 
explicit  results  are  given  in  Ref.  8.  The  basic  idea  is  the  use  of  an  extension  of  the 
calculus  of  variations  called  Pontryagin's  Maximum  Principle  to  determine  the  optimum 
modulation  under  various  constraints  such  as  the  signal's  peak  amplitude,  total  energy 
and  bandwidth. 

The  equivalent  linear  model  suggests  a  straightforward  yet  versatile  approach 
to  computer  investigations.  This  is  potentially  of  great  value  because  of  the  difficulties 
associated  with  closed  form  results.  In  the  single  target  case  the  computations  needed 
to  obtain  the  covariance  matrix  of  the  estimate  errors  are  given  by  Eqs.  (3.  8),  (3.  9)  and 


* 


This  is  the  problem  which  originally  motivated  our  investigations. 


(5.3)  combined  with  the  appropriate  equivalent  linear  model.  <^(0  and  ^(t).  One  of 
several  possible  computational  flows  is  given  in  Fig.  1.  For  the  multiple  target  case, 
a  possible  computational  flow  for  the  equivalent  model  of  Eqs.  (6.  25)  and  (6.  26),  is 
indicated  in  Fig.  2.  Either  figure  could  be  programmed  in  a  straightforward  manner 
on  a  general  purpose  digital  computer  although  a  hybrid  computer  which  combines 
integration  capabilities  with  the  general  purpose  arithmetic  operations  is  even  more 
attractive.  A  versatile  mechanization  suitable  for  a  wide  range  of  experimental  signal 
designs  is  possible  as  changes  in  the  equivalent  linear  model  are  easy  and  do  not 
affect  the  main  computational  flow. 

The  multiple  target  analysis  gives  the  estimate  accuracy  for  a  particular  set  of 
target  separations.  This  means  some  range  of  parameter  values  must  be  investigated 
in  order  to  obtain  a  true  picture  of  resolution  ability.  The  value  of  a  convenient 
computational  algorithm  for  such  investigations  is  obvious. 

Insight  into  the  basic  effects  of  factors  such  as  modulation,  target  separation  and 
bandwidth  is  a  valuable  quantity.  Although  the  best  basis  for  such  insight  would  be  simple 
closed  form  expressions  for  the  covariance  matrix,  the  author  has  found  the  combina¬ 
tion  of  the  information  matrix  and  equivalent  linear  model  to  be  very  useful.  It  is  hoped 
this  is  not  an  isolated  phenomena. 
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